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, We study decoherence due to random telegraph and 1// noise in Josephson qubits. We illustrate 

differences between gaussian and non gaussian effects at different working points and for different 
protocols. Features of the intrinsically non-gaussian and non-Markovian low-frequency noise may 
explain the rich physics observed in the spectroscopy and the dynamics of charge based devices. 
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I. INTRODUCTION 



Considerable progress has been recently achieved in implementing two-state systems using superconducting nanocir- 
cuits. Time-resolved coherent oscillations have been measured in Josephson qubits[l|, |2j, |3| and signatures of the en- 
tanglement of coupled charge qubits have been observed^. Limitations in the performances arise from various noise 
sources 0, IE IE IE IE iioj - often material and device dependent. Recently evidence of noise due to individual impu- 
rities behaving as bistable fluctuators (BF) has been observed both in spectroscopy and in time resolved dynamical 
evolution pd|. Sets of BFs are responsible for low frequency 1// noise fl3l| . which is the dominant dephasing mecha- 
nism common (to different extent) to all solid state implementations. A variety of behaviors, ranging from broadening 
due to a slow environment 0,0 to relaxation limited decoherence [sj has been also observed, this phenomenology being 
I ' strongly dependent on the particular device and on details of the protocol|2l Il4l Il5| . Explanation of this rich physics is 
beyond phenomenological theories describing the environment as a suitable set of harmonic oscillators[EIEE3- Since 
the physical sources of noise are discrete in nature, attention has been devoted to environments made of collections of 
BFsJEIlll In this work we will compare effects of gaussian and non-gaussian noise, and show results indicating that 
i i ■ discrete noise models potentially explain the experimental features, the main open question being a characterization 
beyond phenomenology of the physics of the noise sources |8|- 

Let us consider the Hamiltonian H = Hq — %£a z , representing a qubit (Hq = — |fl • a) anisotropically^ coupled 
to a noise source, described by £. Sensitivity to noise can be modulated by tuning the operating point, i.e. the angle 
9 between z and fi, the qubit splitting ft being also tunable. For classical noise £ = is a stochastic process, 
whereas for quantum noise £ is an operator of the environment. Decoherence results from the sensitivity of the qubit 
C~"~- ■ to the environment. For instance if coupling is weak, relaxation and dephasing rates are Tj -1 = sin 2 9 S(Q)/2 and 
T 2 X = (2Ti) _1 +T%~ 1 , where T' 2 _1 = cos 2 0S(O)/2 is the adiabatic rate, responsible for secular broadening [lf|. 
The power spectrum S(u>) = f 2 (££)^ appears, therefore the qubit is able to "measure" statistical properties of the 
environment at the level of two point correlations. In this regime the qubit is not sensitive to other details, for instance 
whether a certain Lorentzian line shape is due to a bistable fluctuator (BF) giving rise to Random Telegraph Noise 
(RTN) or to a continuous Gaussian process with the same Lorentzian Noise (GLN) spectrum (see Fig. QJa). In solid 
state devices T 2 _1 may be very large thus invalidating the weak coupling theory. This is the case of a sufficiently 
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slow environment, no matter how small is v/Q. Consider for instance the Lorentzian spectrum S(u>) = -g- 

Rates in Fig. ^b show that T 2 _1 ~ v 2 /7, diverges for 7^0. The technical reason of this failure 17] is that the 
, weak coupling expansion parameter is g oc v/j. The high-frequency cutoff sets the characteristic time scale 7 _1 
of the noise, thus problems are encountered for slow noise, g > 1, since the qubit becomes sensitive to details of 
.£h ! the dynamics of the environment. For instance it will distinguish RTN and GLN having the same S(u>). Indeed 
suppose we average the signal from an ensemble of experiments identical, apart for the uncontrolled preparation of 
the environment. The effect of RTN due to a single BF is to determine two angular frequencies for the qubit, f2 and 
fl' = Q[(v/Q + cos#) 2 + sin 2 ^] 1 / 2 , and the signal will show beats at the frequency fl' — Q. If the BF is very slow, 
switches between the two frequencies produce decay with T1.2 ~ 7 _1 . In this re gim e we can identify g = (0/ — fl) Ay 
from the condition[T8j for beats to be observable, g > 1. Instead for slow GLN 10], a standard model in NMRflq]. 
decay is determined by uncertainty in the preparation of the environment, an effect analogous to the "rigid lattice 
line breadth" 0] . By averaging the phase over a static distribution of the effective bias e — v£o , we find the decay 
of the qubit coherences T(t) = -bx\p (t)| = l n (e J *^ 2 + 2 ° £o cose+e o )1/2 ). The distribution of e Q is gaussian for GLN, 
the standard deviation being a 2 = J ^ S(u>) = v 2 /4. The decay depends on 9 and it is slower for 9 = tt/2 where 
for v <§; fl a simple integration[n| yields T(t) = —\ In (l + ia^t / '(2f2)) . For 9 = the full dynamic problem can be 




FIG. 1: (a) Upper panel: a sample of £(t) for RTN and for GLN with the same S(lj). Lower panel (a y ) for a BF with v/£l — 0.2, 
7/fi = 0.02, has almost the same effect as GLN at 9 w n/2. (b) Relaxation and dephasing times by weak coupling theory (solid 
lines) and behavior for slow RTN and GLN (dashed); (c) Same set as in (a) for 8 — 5n/12 with RTN and GLN. In the inset 
the spectral lines for RTN, showing two peaks and for GLN. 
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FIG. 2: (a) Calculated relaxation (thin solid lines) and dephasing (thick solid) times versus Q/7 at 8 — n/4 for decoherence 
due to an impurity with w/fi = 0.2, Sp — 0.08. Grey lines are the weak coupling result, (b) Bottom: The crossover condition 
7c(#). Top: (a y ) for the same BF v/Q = 0.2, 7/fi = 0.05, tuned from the weak (8 = n/2, gray) to the strong (8 = 7r/4, black) 
coupling regime: the appearance of beats can be modulated with the external bias. 



solved and T(t) = (v/2 r y)( , yt — 1 + e~ 7 *) » (vt) 2 /8. In both cases the signal amplitude is strongly reduced already at 
times t ~ v^ 1 <C 7 _1 giving an apparent decay time independent on 7 and finite. Differences between slow RTN 
and GLN are also apparent in the line-shapes (see Fig.njc). 



II. QUBIT DYNAMICS IN THE PRESENCE OF IMPURITIES 



Decoherence due to discrete noise sources has been studied in Ref.[6|,l3 where the ffuctuator has been modeled by 
a quantum impurity described by the Fano Anderson model. Rather than study higher orders of the perturbative 
expansion in v, a different strategy has been used, namely one may write a master equation (ME) for a four-level 
system system including the impurity level, which is traced out at the end of the calculation. Thus fast components of 
the environment are effectively treated at the ME level, whereas slow fluctuations are accounted for nonperturbatively. 
In practice one has at most to diagonalize a simple 8x8 Lindblad map which implies that the qubit dynamics is 
characterized by combining exponentials showing two relaxation times T\± and two dephasing times T 2 ± (Fig.[2Ja). If 
jCl the dominant rates for t ^> 7 _1 coincide with the weak coupling ones (and the classical GLN ones) and describe 
homogeneous broadening. Instead for impurities with j » 1, both are important since dephasing reflects the 
bistable nature of the slow environment. Relaxation is not so sensitive: for v <C f2, the dominant rate is the weak 
coupling one the other rate being T^ 1 » 7. Rates have the form T 2 ~^_ 1 + i5Q/2 = 7(1 ± a)/2 and T^ 1 = 7(1 ± a r )/2, 
dfl being the splitting between the two spectroscopic peaks induced by the bistable impurity 7]. They can be found 
analytically in the full adiabatic regime 7 < O, which includes the crossover region^- In the stochastic limit for the 

BF dynamics a — [(Sp — ig) 2 + cos 4 (<5#)(l — Sp )] and a r = [l — sin 2 (<56')(l — Sp )] , where S9 is the induced 
angular splitting on the Bloch sphere. Features of the nature of the BF appear if Sfl > 7, i.e. > 1, this criterion 
for "strong coupling" depending on both u/7 and the operating point (see Fig- 13b) ■ Strongly coupled impurities are 
non-gaussian, and determine non-exponential decay. Moreover their effect depends on their preparation, so they are 
non-stationary_and non-Mar kovian, determining memory effects. As a consequence decoherence depends on details of 
the protocols [l5j, a critical feature for 1/ / noise. 
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FIG. 3: Results of stochastic Schrodinger equation, with up to Nbf = 1500 BFs producing 1/f noise. Noise level is appropriate 
to charge devices, S(uS) = IQ-kAEq/uj and A ~ 10~ 6 and it is obtained with rid = 250 BFs per noise decade, with v — 0.02fi. 
Operating point is 8 — tv/2 and f2 = 10 10 Hz. In all the examples we show, relaxation coincides with the weak coupling result 
(thick solid line and dots), (a) Adiabatic noise 7m = 10 Hz: the feedback procedure substantially decreases the transverse 
r^(t) (gray curves), whereas inhomogeneous broadening is well described by the rigid lattice breadth theory (dashed lines); 
asymptotic behavior from Makhlin-Shnirman theory is also reported (thin line), (b) Nonadiabatic noise 7a/ = 10 11 Hz: the 
main effect is decoherence during time evolution, determining T2 w Tj. 

III. 1/F NOISE AT OPTIMAL WORKING POINT 

Ideally quantum protocols assume we can measure individual members of an ensemble of identical (meaning that 
preparation is controlled) time evolutions of the qubit. The environment limits control of the preparation. At most 
one can recalibrate, for each individual member, the collective variable £. With this feedback scheme error would be 
due solely to decoherence during time evolution. In actual experiments lack of control on the environment determines 
additional defocusing of the signal, analogous to inhomo gene ous broadening in NMR0, which depends on the 
statistics of the environment at appropriate low- frequencies [l5j . 

We model 1/f noise with a set of impurities switching at rates 7*. For charge based devices they are coupled with 
the qubit via the total polarizing charge, v£ = ViUi, where rii — 0, 1. The standard assumption[l2| of a distribution 
of 7j with P(-f) oc I/7 for 7 G [7m,7A/]|I2 leads to the spectrum S(oj) = Y,i \ v l (1 ~ ^v ) whicn is V/ 

at frequencies 2-nf G [lm, 7m]- In other words, impurities with both large and small g are present. Decoherence 
for this model has been studied for 9 — 0, where exact solutions are available even for non- linear environments |(|. 
The feedback protocol results not to be sensitive to very slow impurities, having 7, < 7* ~ v/10, which sets an 
effective intrinsic low- frequency cutoff 7*. This effect is due to the nongaussianity of the environment, and to the 
fact that BFs with gi ~ 1 determine strong decay. Instead BFs with gi 1 behave as an environment of quantum 
harmonic oscillators, being sensitive only to the amplitude of noise oc n ( iv 2 . In absence of recalibration, inhomogeneous 
broadening corresponds to a proper averaging of the initial conditions of the BFs and one may prove that 7* moves to 
7* ~ min{iT/10, i^ 1 }, where t m is the overall measurement time of the experiment and a large number of repetitions 
is assumed. 

The most effective implemented strategy^ for defeating 1/f is to tune optimally the working point, which in our 
case means setting 9 = n/2. Here the qubit splitting is less sensitive to bias fluctuations, or equivalently the lowest 
order adiabatic rate T' 2 X vanishes. In other words, operating with 9 modifies the parameters gi and the quantitative 
characterization of the spectrum in terms of 7*, the main effect being that smaller gi means smaller effectiveness of 
noise. At first sight one expects that differences between a BF-1// environment and a set of oscillators with 1/f 
spectrum (G-l//) decrease by operating at 9 = tt/2, although they do not necessarily disappear^. Decoherence 
at the optimal point for a G-l// environment has been recently studied by Makhlin and ShnirmanjjJ, combining 
the adiabatic approximation and diagrammatic perturbation theory. Instead we study a BF-1// environment by 
solving numerically the stochastic Schrodinger equation. We checked that very slow BFs are ineffective, so we choose 
7 m = 10 5 Hz. We consider first an adiabatic environment, jm <C ^ (Fig- El a ) studying relaxation via (a x (t)). 
The resulting rate reproduces the weak coupling result, in both protocols with and without recalibration of the 
environment. Instead dephasing is much faster and recalibration is able to reduce defocusing effects. Inhomogeneous 
broadening can be estimated with the rigid lattice line breadth formula. If one assumes Np l to be large enough that 
the initial £0 = J^S 1, Ujn,(0) is gaussian distributed, one finds cr e 2 = v 2 Nbf /4 = 16EcA\n ('yAi/jm), where physically 
7 m = 7* = l/t m . This simple result accounts for the initial reduction of the signal amplitude, and coincides with 
the short-time behavior of the diagrammatic theory^. Adding faster BFs to the environment increases relaxation 
more than dephasing. In our example (Fig.|3]b) decoherence is limited by relaxation T2 ~ T\, this latter being due 
to the fast part of the spectrum u) ~ Q. Inhomogeneous broadening does not reduce drastically the amplitude and, 
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FIG. 4: (a) (cr y ) at 6> = tt/2, SI = 10 10 Hz. The effect of weak adiabatic 1/f noise (light gray line) (7 G [10 5 , 10 9 ] Hz, uniform 
v — 0.002 fi, rid = 250) is strongly enhanced by adding a single slow (7/fi = 0.01) more strongly coupled {vq/Q. = 0.2) impurity 
(black line), which alone would give rise to beats (gray and dashed line), (b) The latter two cases display a characteristic 
behavior of the Fourier transform, shown for the single impurity alone (dashed line) and for 1/f noise plus impurity (solid line). 



as expected, dephasing is underestimated by adiabatic approaches 

Finally, we show that even a single impurity on a 1/f background may cause a substantial reduction of the 
amplitude. This strongly poses the problem of reliability of charge based devices. Effects of realistic BF distributions 
were pointed out by Galperin et al.[2j- Using reasonable parameters, we obtain that indeed an additional BF has a 
substantial effect (see Fig. QJ, whose signature is an asymmetric Fourier transform double-peak, which is similar to 
recent spectroscopy observations fi"l) . 
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